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Abstract: We identify a region of F (R) = gravity without external sources which is isometric 
to the spacetime of colliding plane waves (CPW). From the derived curvature sources, A*' (A'' > 1) 
measures the strength (i.e. the charge) of the source. The analogy renders construction and collision 
of plane waves in F (R) = i?^ gravity possible, as in the Einstein-Maxwell (EM) theory, simply 
because R = 0. A plane wave in this type of gravity is equivalent to a Weyl curvature plus an 
electromagnetic energy- momentum- like term (i.e. 'source without source'). For = 1 we recover 
naturally the plane waves (and their collision) in Einstein's theory. Our aim is to find the effect of an 
expanding universe by virtue of F{R) = R^ on the colliding gravitational plane waves of Einstein. 

o 

F (R) gravity in which the Einstein-Hillbert action (i.e. F (R) = R) gets corrections by an arbitrary function of 

■ scalar curvature R has become mucb popular in recent years. Such an interest stems from cosmological observations 

■ and attempts to interpret the accumulated data. The accelerated expansion of our universe with dark energy / matter, 
phantom fields, quintessence and cosmological constant can all be addressed under a common modified theory entitled 
as F (R) gravity. Dependence on the single (and simplest) Ricci scalar R becomes an advantage when compared with 
theories depending on the higher order invariants. The F (R) gravity, even in the absence of external sources is highly 

I ■ non-linear, but recent studies show that many exact solutions have been obtained [l| (for the review of the subject), 
^jj^ No doubt, the freedom of adding suitable energy momentum in addition to those arising from curvature makes the 
I— I. theory further complicated, but interesting. 

In this paper we wish to draw attention to another aspect of F (i?) gravity which is local equivalence (isometry) 
. with the spacetime of colliding plane waves (CPW) (for the review of the subject see 0) . It was shown first by 
^ • Chandrasekhar and Xanthopoulos (CX) in 1986 Q that part of black hole spacetime is locally isometric with the 
[ spacetime of CPW [3] . 

■ For black hole spacetimes with single (event) horizon, such as Schwarzschild, the region of isometry is inside the 
I horizon where the geometry admits two space-like Killing vectors. For black holes that admit double horizons (i.e. 
. r = r_ = the inner and r = r+ = the outer horizon) such as Reissner-Nordstrom (RN) and Kerr, the dynamical 
" isometry region occurs in between the two horizons, r_ < r < r+. This property has been used in the past to generate 

y—i . CPW solutions from the known solutions of black holes 4] . Specifically, this technique has aided in generating regular 
^-H ■ classes of CPW solutions due to the fact that the singularity of the black hole is removed in the isometry process. 
^ ] Let us note that the intrinsic connection between black holes and plane waves is not limited by the observation in . 
• ■ Much earlier Penrose proved a Theorem Q in this connection: " Every spacetime admits a plane wave spacetime as a 
[ limit" , and naturally black hole spacetimes are not exempted from this statement. For a class of extremal black holes 
we have shown even further, that the near horizon geometry coincides with that of CPW '6^. Being motivated by these 
past studies it is our aim to point out that similar isometrics exist also in F (R) — R^ gravity and in a larger class of 
F (R) theories which admit RN solution of Einstein's theory as a particular solution. To be on the safe side, however, 
in the satisfaction of F (R) equations in the limit i? — )■ 0, we restrict ourselves to the F (R) = i?^ case. In this study 
we consider a 4 -dimensional spherically symmetric F (R) gravity solution without external sources. That is, the 
curvature creates its own source through the intrinsic nonlinearity. This derived source, or more popularly, "source 
without source" , turns out to be an electromagnetic-like source and the geometry is defined modulo the RN geometry. 
Under a suitable transformation this part transforms into the spacetime of CPW. The set of boundary (junction) 
conditions for the CPW is adopted from the Einstein-Maxwell (EM) theory in which the appropriate conditions are 
those of O'Brien and Synge These guarantee that non-physical terms, such as the product of two Dirac delta 
functions, or any other undesired combinations do not arise at all. It should be added that a general analysis of 
junction conditions in more general F (R) gravity other than F (R) = R^ , deserves a separate study. 



•Electronic address; lhabib.mazhari@emu.edu.tr | 
tElectronic address: mustafa.halilsoy@emu.edu. tjj 
Electronic address: ,tayabeh.tahamtan@emu.edu.tr| 



2 



By construction our solution of CPW admits a horizon following the collision, in place of a spacetime singularity. 
Whether the horizon is stable or not, and if stable / unstable, against what type of perturbations, remain open to 
which we shall not address in this paper. Further, by going from the interaction region (Region IV) to the incoming 
regions (Regions II and III) we identity, in analogy with taking the Penrose-limit [5] , the incoming plane waves prior 
to the instant of collision. Region I is flat or the no-wave zone. Extension into the incoming regions give us a class 
of plane waves in F (R) — gravity which consist of a gravitational (Weyl) and an energy-momentum part derived 
from the curvature. In this picture {N — 1) plays the role of the electromagnetic charge so that we assume iV > 1. 

Similar analysis can be considered for different versions of F (R) gravity, provided the function F (i?) is known 
explicitly. Let us note that in most of the cases considered in cosmology a closed form(i.e. analytical) representation 
for F (R) is not available. From physics standpoint a plane wave in general relativity is a representation of an infinitely 
boosted particle under proper limiting procedure 8] . At such high speeds a particle loses its matter property and turns 
into a light-like wave. For this reason in the large hadron collision (LHC) experiments the particle can alternatively and 
approximately be considered as a wave. From this token it should be added that the almost plane wave representation 
of an infinitely boosted proton, due to its internal structure (i.e. colors, gluons, quarks, etc.) is yet to be found. Our 
boosting process covers so far only point-like objects without internal structure. In this regard what is available at 
our disposal is no more than a coarse approximation to the exact case. Another important item is that the horizons 
formed in CPW are unstable against arbitrary perturbations @. This may be the classical justification that micro 
black holes formed at high (TeV scale) energies are not stable / detectable. From the quantum picture the probability 
of black hole (or any other objects, such as mini wormholes) formation is non zero. However, if these objects are to 
be detected classically their stability must be verified. In this case the stability arguments of horizons in CPW may 
offer some hints. It is believed that modified theory of gravity such as F (R) with its rich content will give tangible 
results. 

Organization of the paper is as follows. Section II summarizes F (R) gravity with its equivalence to the space of 
CPW. Section III discusses the properties of our colliding wave spacetime. Comparison of Einstein and F{R) = R^ 
theories takes place in Section IV. Conclusion and Discussion in Section V completes the paper. 

II. F (R) = GRAVITY AND LOCAL EQUIVALENCE WITH THE SPACE OF CPW. 

We consider the action in the form (SttG = 1) 



^Jd^xV^[F{R) + L^], (1) 



where the function F (R) will be determined by the field equations in general and L„i stands for the matter Lagrangian. 
In this paper we shall confine ourselves to the particular case in which Lm = 0, so that any 'source', T^^ will be 
determined by the curvature. Further, we shall follow an easier route to assume F (R) — R^ a priori, and use it to 
find the curvature sources. With these assumptions Einstein's equations are still valid as 



F'{R) 



i^M. (F (i?) - RF' (R)) + (V^V, - g^,n) F' (R) 



(2) 



where F' (R) = |^ and □ — V^V^ . When the geometrical components of G^i^ are substituted into the foregoing 
equation it yields 



Rf^.F' (R) - V^V.F' (i?) + (of' (R) - ^F {R}^ g^, = 0. 



(3) 



By taking the trace of the field equation one obtains 

F' (i?) R + iUF' (R) - 2F {R) = 0. (4) 
Now, assuming a spherically symmetric metric ansatz 

ds^ = f (r) df - — - - (de^ + sin^ 6 dip^) , (5) 
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in which / (r) is the only function to be found and the choice F (i?) = yields [T^ 

lr/" + 2/'\ 1 



N 



r 



N 



rf' + {f-l) \ Ip 

7^ 



= 0, 
= 0, 



i?^ ( 1 



N 



dr 



(6) 
(7) 
(8) 



The latter equation admits cither = 2 or i? = which is the case of our interest here. The other equations also 
confine the parameter to > 1 in order to be satisfied. Therefore setting R — admits a solution for / (r) as 



„ , , 2m a 



(9) 



Here m and a are two integration constants. Throughout the analysis we shall make the choice <j ^ 0, since otherwise 
/ (r) reduces to that of Schwarzschild, and also our choice will be that of TV > 1. This form of / (r) is easily identified 
as the RN black hole where m and a are the mass and square of the charge, respectively. In the following analysis 
we shall make the simple choice m = 1 without loss of generality. This class of solutions can be considered in three 
difi^erent cases, depending on the horizon(s), namely 



This admits the roots 



for 



2r±+a = 0. 



r± = 1 ± a, 



(10) 



(11) 



(12) 



We can identity three different classes: 

i) a = 1, or (7 = 0; This brings us back to Einstein's gravity, so it is not interesting for our purpose. 

ii) a = 0, or (T = 1: This choice corresponds to the extremal RN whose near horizon geometry is the Bertotti- 
Robinson (BR) spacetime. In other words, we have 



Now, making the transformation 



f{r)=[l-- 
1 r 



1 + A7~, 



t 

A' 



in the limit, A — ?► 0, with the successive replacements f — and t = T yields for the line element 



1 

R^ 



{dT^ ~dR^)- de^ ~ sin^ B difi^ 



(13) 



(14) 



(15) 



This is the well-known BR spacetime of EM theory which can be transformed easily into the space of CPW by 
a well-established transformation The resulting spacetime will be the Bell-Szekeres space time for colliding 
electromagnetic shock waves [ll|. For this reason we shall not be interested in this particular case of a = 0, either. 

iii) < a < 1 (or < cr < 1): In this paper we shall be interested in this particular class since (i) and (ii) will 
not add much to our interest. From now on we can make the choice < cr < 1 (0 < a < 1), which is the consequence 
of choosing N > 1 from Eq.s (6-8) . Our line element (5) reads now 



= f 1 - - + ^ ] df2 _ 



sin^ 9 dip^ 



(16) 
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In between the horizons 

l-a = r_ < r < r+ ^ 1 + a, (17) 
t becomes space-Hke and r becomes time-hke. Now let us apply the transformation 

r = 1 + ar, cos6 = a, <f = y, t = x, (18) 

to the new coordinates (T,a,x,y) from the old ones {t,r,9,Lp). Next, we transform once more from (t, cr) to the 
null-coordinates (u, v) by 

T — sin {a u + b v) , a — sin {a u — b v) , (19) 

in which (a, b) are constants that we introduce for convenience. These measure the magnitudes of incoming energies 
and in case that we choose a = 1 (6 = 1) we normalize those strengths to unity. After an overall scaling of coordinates 
X, y and ds^ by the constant 2ab the line element takes the form 

ds^ = (1 + a rf {2du dv - Sdy"^) ^dx^, (20) 

(1 + a r) 

in which we have abbreviated A — 1 — and S = 1 — a^. So far, we have only applied coordinate transformations on 
our original metric (5) until we reach at (20). It is a remarkable observation at this stage that the metric (20) can be 
interpreted as a CPW metric under the substitutions 

u->-u0{u), v->-v0{v), (21) 

where 9 (u) and 9 (v) are the Hcavisidc step functions corresponding to the two different null boundaries u = and 
V ~ 0. The adopted O'Brien-Syngc [7] junction conditions guarantee in this process that no extra sources arise at the 
null boundaries. Otherwise the metric can not be interpreted as a CPW spacetime and the idea should be discarded. 

III. PROPERTIES OF THE COLLIDING WAVE SPACETIME 

Upon insertion of the step functions our CPW reads 

ds^ — (1 + a sin [au 9 (u ) + bv 9 {v)]f {2du dv - cos^ [au 9 (u) ~ bv 9 {v)] dy^) - 

cos^ [au 9 {u) + bv9 jv)] ^^^^ 
{1 + a sin [au 9 (u) + bv 9 (w)])^ 

With the choice of null-basis 1— forms of Newman and Penrose (NP) (in which the step functions present) 

I = (1 + ar) du, 

n = (l + ar)dw, (23) 

[6 ,^ . , . [a dx 
m = ^-il + ar)dy + ^^-J^^^ 

all non-vanishing spin-coefficients and NP quantities are tabulated in Appendix A. Fig. 1 depicts the details of the 
space time of CPW. The Region I (u < 0, w < 0) is flat. The Region II (u > 0, w < 0)has the wave profile 



(24) 



5 



which is the superposition of an impulsive and a shock gravitational wave. Since our class of solutions is valid 
for < a < 1 the impulse wave has trailing shock waves and both components of impulsive and shock waves are 
indispensable. The only surviving energy-momentum component in Region II is 

<^2^{u)^—^-^^—-L. (25) 
(1 + a sm (a u)) 

The incoming Region III has the same structure as Region II with a — > 6, u — > w so that the non- vanishing components 
are '^q (u) and $00 (i') (see Appendix A3). After collision, in the interaction region (Region IV) (m > 0, u > 0) we 
have the additional components of ^2 (w;U) and $02 (u,w) = <I'20 (u,f) which didn't exist in the incoming regions. 

It is observed from the NP invariants that the interaction region (u > 0, t; > 0) has no singularities (i.e. A (u, w) — 
(1 + a sin (au + t'w )) 7^ 0. On the null-boundaries, however, (i.e. au — ^, u = and hv = -j, m = 0) there are null- 
singularities since the Weyl components diverge there. In the interaction region (it > 0, w > 0) none of the curvature 
invariants diverge, the metric function g^x , however, vanishes which indicates a coordinate singularity. It can be 
recalled that before applying the coordinate transformation g^x was gu whose zero gave rise to the horizon. The 
hypersurface sin(au + bv ) = 1, is the horizon that forms as a result of colliding waves. From the CPW experience 
in the EM theory it is known that such horizons are unstable against perturbations That is, addition of any small 
physical field gives rise to a back reaction effect which makes the horizon to transform into a true spacetime singularity. 
By extension beyond the horizon (i.e. au + bv > ^) it is evident, for example somewhere within ^ < au + bv < 
that we face a curvature singularity. It is expected that the same behavior takes place also in F (R) gravity. Our 
example of solution is simple enough to simulate the CPW in EM theory. Any solution that will be found in this 
regard must inherit a more violent wave interaction and stronger singularity. In the interaction region (m > 0, ?; > 0) 
we observe that off-thc null boundaries the Dirac delta functions vanish and only the shock terms survive. Under this 
condition the Weyl curvature components satisfy 



9*^ = 1'o*4, (26) 

which implies that the spacetime is of type-D, expectedly since after all it is a transform of a type-D spacetime. In this 
spacetime it can easily be shown that the massive scalar wave equation and Hamilton- Jacobi equations can easily be 
solved by separation of variables. Let us add that the null coordinates (u, v) is not an appropriate coordinate system 
for this purpose. A more appropriate coordinate system is given by (r, a, x, y), so that the scalar field equation is 

(□ + m) $ (r, (7, a;, y) = 0, (27) 
where □ is the covariant Laplacian and m is the mass of the particle. We seek a separability in the form 



$ (r, a, x,y) = T (r) S (a) e*^""+™«^ (28) 

where /ip and vq are constants and the functions T (t) and S (cr) are functions of their arguments. The decoupled 
equations are reduced into quadratures 



(AT^r + (^^^ -ko)T = 0. (29) 

with the separation constant ko , A (r) is defined in the Appendix A2) and sub (r) / (cr) indices imply partial 
derivatives. In the next Section we shall make a comparison between general relativity and F{R) gravity results in 
both the NP quantities and the restricted geodesies. 



IV. COMPARISON BETWEEN F{R) = R AND F{R) ^ R'^ (N > 1) 



A. The Ricci and Weyl tetrad scalars 



All the related tetrad scalars are given in the Appendix (A3-A5). For a = 1 they reduce to F{R) = R which are 
the results of general relativity. Now, in order to compare with the case of F{R) — i?^ we set 



a = 1 — e 



(30) 
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in all quantities to the first order with e > and ~ 0. Upon substitution and expansion in first order of e we obtain 
the following forms 

= + + O {e') (31) 

^r, = + + o (e2^ (32) 

in which <^^^^ and ^f'^ are the Ricci and Weyl components in the Einsteins general relativity. The ^jj-* and ^'f'^ are 
naturally the first order contributions from the F{R) gravity. As a result the expressions (A3-A5) are obtained as 
follow 



^(0) ^ -a6{u) ^ 3a\^,(o)_ -b5{v) ^ 3b^ ^ 

^ cos6w (1 + sin6z;)'^ S ' " cos au (1 + sin att)^ ' 

<> = = = = ""fr/.-^' . (33) 

^ ^ ^ cosow (1 + smot;) 

^(1) _ -a^(u)(3sin6t;- 1) _^ 3a2(3sinii'-2) 

^ cos foil (1 + sin fef)"^ 5]^ 

(1) — &(5 (f) (3 sin ait — 1) 3a^ (3 sin — 2) 
cosaw (i + smau) ^ 

in which S = 1 + sin K and K is defined in Appendix A. We have also 

Zm („) = ^ZL™ ^ zW {u) = ^( , (34) 
COS au cos au \ sin au + 1 / 

Zio) (,) = ^Zlfi^l^, ^(1) (z;) = ^ f ^^"^^-^ 
cos bv cos 6w \ sin bv + 1 

By choosing u < 0/w < in the above expressions we find the differences between the incoming plane waves between 
the two theories. This may be interpreted as the distortion of Einstein's plane gravitational waves in the expanding 
universe model due to the F{R) = R'^ gravity. These expressions complete the amplitudes of F{R) = R'^ gravity in 
comparison with the general relativity at the first order. 



B. Deviation of Geodesies 

We consider the geodesies projected into the {u, v) plane (i.e. x =cons., y =const.) and for simplicity we fix the 
constants a = b = 1. The system can be described by the Lagrangian 

L{u,v) = {1 + asm.{u + v))\/v' (35) 

in which v' = 3^ , so that v is parametrized in terms of u. With this much simplification the geodesies equation follows 

as 

^// ^ 2a^^ (1 - 2v') cos [u + v) 
1 + a sin (u + v) 

Obviously an exact solution follows for the choice v — ^u, but in general the differential equation is highly nonlinear 
which can be expressed in terms of elliptic functions. Since our aim is to compare geodesies with general relativity 
we take again a = 1 — e, as we did above and make expansions in terms of e. We obtain to the first order in e the 
following differential equation 

v" = 2v' (1 - 2v') ( 1 - --4 -) . (37) 

^ ^ 1 + sm (u + w) V l + sm{u + v)j ^ ' 

Fig. (2) displays the numerical plots of Einstein's (e = 0) and geodesies for the small parameters e = 0.2. The initial 

speeds v' = 3^ is varied accordingly. The geodesies take place within the region u + v< ^, restricted by the horizon 
shown by the broken line in Fig. (2). 
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V. CONCLUSION AND DISCUSSION 



The pathological problems associated with the singularities and horizons in colliding waves were exclusively discussed 
in the past (see 0] and references cited therein) so that we didn't feel the necessity of repeating them here. Our aim 
is to draw attention to a well-known local isometry between CPW and black holes that exist also in the nowadays 
fashionable F {R) gravity. We have chosen the particular case of F (R) = R^ , (TV > 1) as an example in which the 
scalar curvature R — 0. Our solution is single parametric with parameter a=^/l — a . The employed solution is locally 
isometric to a portion of the Reissner- Nordstrom (RN) geometry and effectively our F (R) = i?^ gravity becomes 
Einstein-Maxwell (EM)-like. The physical interpretation, however, is entirely different. There are no separate Maxwell 
equations for instance. In this way from the local isometry we construct CPW solutions from a particular F {R) gravity. 
By interpolating into the incoming regions we obtain a class of plane waves (i.e. the Penrose limit) in this theory. In 
general, the Weyl curvature of these plane waves are in the form of (Impulsive) -I- (Shock) waves. Its energy-momentum 
in all regions consists of pure shock terms, as it should be. From the EM analogy the boundary (junction) conditions 
that has been adopted are those of O'Brien- Synge which eliminate any spurious boundary currents. We remark 
that for more general (and further nonlinear) structure of F (R) theories new boundary conditions are required. We 
admit that our metric of CPW is the simplest one in which by virtue of i? = all pertinent F (R) equations are 
automatically satisfied. We find all deviations of the Weyl and Ricci scalars from the Einstein's gravity and compare 
also the geodesies to the first order. In analogy with the colliding electromagnetic and Yang-Mills jl^ shock 
waves our interaction region here also forms a horizon through which the spacetime can be extended analytically 
14. As the next step we may consider the case of i? = i?o =constant. This will give rise to a cosmological constant 
1^ whose colliding wave interpretation can not be considered free of null shells p^ . Finally we wish to comment 
on the importance of the collision problem presented in this paper. Infinitely boosted particles can be represented 
by plane waves in general relativity. Head-on collision of such plane waves simulates high energy particle collisions. 
For this reason production and detection of exotic objects such as microscopic-black holes / wormholes becomes a 
reality provided their classical stability analysis is verified. No doubt the richer structure of F (R) gravity and the 
mathematical analysis of horizon stability incorporated in CPW can be helpful in this context. 
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Appendix A 

a) The non-zero spin coefficients are (for < a < 1) 
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where we have used the abbreviations 



a bO (v) cos K , . ^ , 
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a aO (m) cos K 
2^2 



a0 (u) , 

II = -^^ (tanii^ + tanL) 

p = —Y-{ta.nK-ta.nL) 
A = ^A^(Z-tani) 
a = -^(Z + tani) 



A = l + asinK, (A2) 
K = aue{u) + bv e{v), 

L = a u 9 (u) — b V 9 (v) , 
s'mK + a{l+ cos^ K) 



Z = 



AcosK 



b) The non-zero Ricci components are 



^oo{u,v) = (A3) 

, , a^9{u)(l-a^) 

ab0{u)9{v){l-a^) 

^20iu,v) = ^Q2{U,V) = 

c) The non-zero Weyl components are 

^ , , aab9 (u) 9 (v) (a + sin K) 

^2{u,v) = (A4) 

. X b5(v) , , 2,ab'^9{v)(a^s\nK) 
= -^]4^ (^(«) + tan («))) + "^-^^ ^, 

in which (5 (m) and (5 (v) stand for the Dirac delta functions and Z (u) and Z [v) are 

sin (a^i^^ (u)) + a (l + cos^ (au(? (u))) 
^^"-'^ cos(aM0(u))(l + asin(aM0(M))) ' ^ 
Z (w) = Z (m ^ , 
(Z(0) = 2a). 

Figure caption 

Figure 1: The (it, u) plane of CPW consists of four disjoint regions: Region I (u < 0,w < 0) is flat. Regions II 
(w > 0, < 0) and III (u < 0, ii > 0) are the non-flat incoming regions and Region IV (w > 0, > 0) is highly curved 
Interaction region. Prom the right (Region II) gravitational wave (u) (accompanied with the energy components 
$22 (w)) collides with ^'q (w) from the left (superposed with the energy term $00 (''^)) at u = = w. The collision 
develops the Region IV (the Interaction Region u > 0, ii > 0) in which the nonzero curvature (and energy) components 
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are shown. The hypersurface au + bv = ^ is a horizon whose intersections with the miU boundaries (i.e. au ^ ^, u = 
and bv = = shown with dark points) arc null singularities. The lines (i.e. au = ^, ?; < and bv — < 0) are 
horizons of the incoming regions. 

Figure 2: Particular geodesies in the {u,v) plane of colliding waves according to the differential equation (35). 
Curves labelled by A are the geodesies of general relativity while B's are the F{R) geodesies when the infinitesimal 
parameter is taken as e = 0.2. Depending on the initial speed v' = the difference of two classes is revealing. In 
particular, for v' > 1 the difference becomes significant while approaching the horizon given hy u + v = ^. The only 
geodesic that remains invariant in going from F{R) = R to F{R) = R^ is the case v' = 0.5 which is evident from 
Eq. (36) in the text. In effect, the expansion of the universe due to F{R) gravity manifests itself in the deviation of 
geodesies as displayed here. 
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